Scattering approach to parametric pumping 
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A d.c. current can be pumped through a quantum dot by periodically varying two independent 
parameters Xi and X2, like a gate voltage or magnetic field. We present a formula that relates the 
pumped current to the parametric derivatives of the scattering matrix S(Xi,X2) of the system. 
As an application we compute the statistical distribution of the pumped current in the case of a 
chaotic quantum dot. 

PACS numbers: 72.10.Bg, 73.23.-b, 05.45. +b 
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An electron pump is a device that generates a d.c. cur- 
rent between two electrodes that are kept at the same 
bias. In recent years, electron pumps consisting of small 
semiconductor quantum dots have received, considerable 
experimental and theoretical attentionBilil A quantum 
dot is a small metal or semiconductor island, confined by 
gates, and connected to the outside world via point con- 
tacts. Several different mechanisms have been proposed 
to pump charge through such systems, ranging from a 
low- frequency modulation of ga±&, woltages in combina- 
tion with the Coulomb blockadeclfl'til to photon-assist 
transport at or near a resonance frequency of the dot. I 
Their applicability depends on the characteristic size of 
the system and the operation frequency. 

Most experimental realizations of electron pumps in 
semiconductor quantum dots made use of the princi- 
ple of Coulomb blockade. If the dot is coupled to the 
outside world via tunneling point contacts, the charge 
on the dot is quantized, and (apart from degeneracy 
points) transport is inhibited as a result of the high en- 
ergy cost of adding an extra electron to the dot. Poth- 
ier et al. constructed an electron pump that operates at 
arbitrarily low frequency and with a reversible pump- 
ing directions The pump consists of two weakly coupled 
quantum dots in the Coulomb blockade regime and oper- 
ates via a mechanism that closely resembles a peristaltic 
pump: Charge is pumped through the double dot array 
from the left to the right and electron-by-electron as the 
voltage U\ oc sin(wt) of the left dot reaches its minima 
and maxima before the voltage U2 oc sin(wt — cj>) of the 
right one.13 The pumping direction can be reversed by 
reversing the phase difference (j> of the two gate voltages. 

A similar mechanism was proposed by Spivak, Zhou, 
and Bcal Monod fo* an electron pump consisting of single 
quantum dot onlyD In this case a d.c. current is gener- 
ated by adiabatic variation of two different gate voltages 
that determine the shape of the nanostructure, or any 
other pair of parameters X\ and X2 , like magnetic field 
or Fermi energy, that modify the (quantummechanical) 
properties of the system, see Fig. |l|a. The magnitude 
of the current is proportional to the frequency u> with 
which X\ and X2 are varied and (for small variations) 
to the product of the amplitudes 8X\ and 8X2- The di- 
rection of the current depends on microscopic (quantum) 



properties of the system, and need not be known a priori 
from its macroscopic properties. As in the case of the 
double-dot Coulomb blockade electron pump of Ref. |, 
the direction of the current in the single-dot parametric 
pump of Spivak et al.ta can be reversed by reversing the 
phases of the parameters X\ and X%. An important dif- 
ference between the two mechanisms is that a parametric 
electron pump like the one of Ref. ^ does not require that 
the quantum dot is in the regime of Coulomb blockade; it 
operates if the dot is open, i.e. well coupled to the leads 
by means of ballistic point contacts. Experimentally, an 
electron pump in an op<in quantum dot has been real- 
ized only very recently.tHI A measurement of the pumped 
current provides a promising tool to study properties of 
open mesoscopic systems at zero bias or at zero current. 

In this paper we consider a parametric electron pump 
through an open system in a scattering approach. Our 
main result is a formula for the pumped current in terms 
of the scattering matrix S{Xx,X%). Such a formula is 
the analogue of the Landauer formula, which relates the 
conductance G — SI/8V of a mesoscopic system with two 
contacts to a sum over the (squares of) matrix elements 
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The indices a and (3 are summed over all channels in the 
left and right contacts, respectively, and SV is the applied 
voltage. For the case of the parametric electron pump, 
where two parameters X\ and X2 are varied periodically, 
8X\ (t) — 8X1 sm(uit) and 8X%(t) — 8X2 sin (ojt — </)), we 
find that the d.c. component of the current / depends on 
the derivatives dS a 0/dX, 
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Like the Landauer formula, Eq. (|2|) is valid for a phase 
coherent system at zero temperature and to (bi)linear re- 
sponse in the amplitudes 8X1 and 8X2- [The nonlinear 
response is given by Eq. (^) below.] It captures both a 
classical contribution to the current and the quantum in- 
terference corrections. Quantum corrections can be im- 
portant in the mesoscopic regime, especially if there is 
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FIG. 1. (a): A quantum dot with two parameters X\ and 
X2 that describe a deformation of the shape of the quantum 
dot. As X\ and X2 are varied periodically, a d.c. current I is 
generated, (b) In one period, the parameters X\ (t) and X2 (t) 
follow a closed path in parameter space. The pumped current 
depends on the enclosed area A in (Xi, A2)-parameter space. 



no "classical" mechanism that dominates the pumping 



process.El'£3 Eq. (|^) is valid to first order in the frequency 
lu. This is sufficient if the period r = 2tt/lu is much 
larger than the time particles spend inside the quantum 
dot. For such low frequencies, we can assume that equi- 
librium is maintained throughout the pumping process. 
The scattering matrix formula does not capture effects of 
order uJ 2 (or higher) that rely on the existence of a. non- 
equilibrium distribution inside the quantum dot.Q The 
existence of a scattering approach to parametric pumping 
allows us to borrow from the vast literature dealing with 
scattering matrices of disordered and chaotic microstruc- 
tures and their parameter dependence,t3 and to directly 
relate the pumped current to other transport properties 
like e.g. the conductance. 

The system under consideration is shown schematically 
in Fig. ^a. It consists of a quantum dot, coupled to two 
electron reservoirs by ballistic point contacts. The two 
electron reservoirs are held at the same voltage. Two 
external parameters Xi (t) and X2 (t) of the dot are varied 
periodically, see Fig. |a. They can be e.g. the voltage of a 
plunger gate, parameters that characterize the shape, or 
a magnetic field. The two point contacts, which have N 
channels at the Fermi level Ep , are labeled 1 and 2. The 
scattering matrix S of the system has dimension 2A*" x 2N 
and is a function of the parameters X\ and Xi . Since the 
system is well coupled to the leads, the charge is no longer 
quantized, the Coulomb blockade is lifted, and to a first 
approximation, we can use a picture of non-interacting 
electrons Jifl 

Starting point of our thecoi is a formula due to 
Biittiker, Thomas, and PretreO for the current in the 
contacts 1 and 2 that results from an infinitesimal change 
of a parameter X: For a small and slow harmonic vari- 
ation X(t) = X + bX u e %u3t , the charge 6Q(m) entering 
the cavity through contact m (m — 1,2) reads 
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The index a is summed from 1 to N for contact 1 and 
from N+ 1 to 2N for contact 2. rljhe quantity dn(m)/dX 
is the emissivity into contact m.t-3 Eq. (^) is valid to first 
order in the frequency oj and assumes that the scattering 
properties follow the time-dependent potentials instanta- 
neously After Fourier transformation one obtains 



dX 



(4) 



Similarly, for a simultaneous infinitesimal variation of 
two parameters X\ and X2, the emitted charge 6Q(m,t) 
through contact m is (m = 1, 2) 



SQ(m,t)=e^^-6X 1 (t) 
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Next, we consider a finite variation of both parameters 
Xi and Xi. The total charge emitted through contact 
m in one period r = 2ix jio is found from integration of 
Eq. (^) to X\ and X2, bearing in mind that the scatter- 
ing matrix S and hence the emissivities dn(m) / dX\ and 
dn(m)/dX2 are functions of X\ and X2, 
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In one period, the pair of parameters X\(t) and X^ii) 
follows a closed path in the (Xi , X2) parameter space, see 
Fig. 0b. The total charge expelled from the dot through 
contact m can be rewritten as a surface integral over the 
area A enclosed by the path in parameter space using 
Green's theorem, 
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/ d dn(m) 

dX x dX 2 — K —L 

a \9Xi 



dX 2 



d dn(m) 
~dX 2 dXi 



Note that the surface area A, and hence the transported 
charge, vanish, if the parameters X\ and X 2 vary in 
phase, or with a phase difference n. The surface area 
is maximal if their phases differ by 7r/2. Substitution of 
Eq. (^b|) for the emissivities yields 
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Hence the d.c. current I m through contact m is given by 
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One verifies that I\ = —I2, indicating that no charge 
is accumulated. The response matrices Rx x and Rx 2 
are hermitian 27V x 2N matrices. For the (bi)linear re- 
sponse to the variations of the parameters X\ and X 2 , 
Eq. (^) simplifies to the result (||) quoted in the intro- 
duction. Note that, since the parameters X\ and X2 are 
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dimensionless, the current formula contains no factor h, 
unlike the Landauer formula (|l|). Planck's constant may 
however appear in the typical scales for the parameter 
dependence of the scattering matrix S(Xi, X2). 

Eq. (0) is the main result of this paper. It establishes 
the link between the pumped current / and the paramet- 
ric derivatives of the scattering matrix S. Several quali- 
tative observations can already be reached on the basis of 
Eq. (||). First, for a phase coherent quantum system, the 
out-of-phase variation of any pair of independent param- 
eters will give rise to a d.c. current to order lo. Second, / 
is not quantized, unlike in the case of the electron, pumps 
that operate in the regime of Coulomb blockade.cl Third, 
if the size of the variations 8X\(t) — 6X\ sin(uit) and 
6X2 (t) = 6X2 sin(wi — <f)) is small compared to the char- 
acteristic correlation scales X\ c and X2 C needed to change 
the scattering properties of the sample, we may neglect 
the Xi and X2 dependence of the integrand in Eq. (||), 
and recover the (bi)linear response formula On the 
other hand, for 5Xj ^ Xj C (j = 1,2), the integrand in 
Eq. (Ii) may have multiple sign changes within the inte- 
gration area A, and the typical value of / is proportional 
to {5Xi5X2X\ c X2c sin (j)) 1 / 2 . [Although the typical value 
of the current scales as (sin (j)) 1 / 2 , the ^-dependence of the 
sample-specific current may be quite random.] 

Like the Landauer formula, the scattering matrix for- 
mula (||) describes both a classical contribution to the 
current and the quantummcchanical corrections. Their 
roles are illustrated below in two examples. First, we 
consider a simple pump in a one dimensional wire. The 
wire contains a tunnel barrier at x — and for < x < L 
a region where the electrostatic potential U can be var- 
ied, e.g. by varying the voltage of a nearby gate, see Fig. 
0a. The Schrodinger Equation for this system reads 



k 2 ^>(x) = f-|j+y(x)j ^(x), 
V{x) =~/5{x) + U9[x(L-x)], 



(9) 



where 6(z) = 1 if z > and otherwise. We pump 
electrons through the system by opening and closing the 
tunnel barrier and raising and lowering the potential U 
as indicated in Fig. [|b. This system operates as a classi- 
cal "peristaltic" electron pump. To find the d.c. current 
/, we compute the scattering matrix S and apply the 
scattering matrix formula (pi), 
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The first term is the classical contribution to the pumped 
current. [Note that the local density of states for this one- 
dimensional system is l/(2-7rfc).] The second term is the 
correction due to quantum interference. 

As a second example, we consider the case where elec- 
trons are pumped through a disordered or chaotic quan- 
tum dot. This application is relevant for the experiments 
of Ref. O. The two (dimensionless) parameters X\ and 
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FIG. 2. (a) An electron pump, consisting of a 

one- dimensional wire with a tunnel barrier at x = and 
an adjustable electrostatic potential U for < x < L. (b) 
Charge is pumped through the wire by varying the height 7 
of the tunnel barrier and the potential U in the following or- 
der: (i) 7 — > 00 (close barrier), (ii) U — > SU (raise potential), 
(iii) 7 — > (open barrier), (iv) U — > (lower potential). 

X2 characterize two different deformations of the shape 
of the dot, see Fig. [l]a. Unlike in the previous example, 
where the pumping mechanism was of a mainly classical 
origin, for a chaotic quantum dot, there is no "classi- 
cal" contribution to the pumped current. The current 
results from quantum interference and its size and direc- 
tion depend on microscopic details of the system. Pump- 
ing occurs because the wave functions near the two point 
contacts are different and strongly parameter dependent, 
so that different amounts of current flow through the two 
contacts if the parameters X\ and X2 are varied. 

For a disordered or chaotic quantum dot, the sta- 
tistical distribution of the scattering matrix S(X\,X2) 
and its dependence [On X\ and X 2 are given by Ran- 
dom Matrix Theory.E3 Within Random Matrix Theory, 
the parameter-dependent Hamiltonian H(X±, X2) of the 
quantum dot is replaced by a large M x M hermitian 
matrix H(X\, X2), 

H(X U X 2 ) = H + M- 1/2 X 1 H 1 +M- 1/2 X 2 H 2 , (11) 

where the matrix elements of H, Hi, and H2 are inde- 
pendently and identically distributed Gaussian random 
numbers. A distinction is made between the cases that 
time-reversal symmetry (TRS) is present [H(Xi, X 2 ) is 
real] or absent [H(Xi, X2) is complex]. We now compute 
the distribution of the pumped current / in the regirpe 
of small parametric variations SXj <C Xj C (m = l,2).li3 
Using the distribution of the matrices Rx ± and Rx 2 from 
Ref. [17], we find that for many-channel leads (N 3> 1), 
the current / is Gaussian distributed with (I) = and 
r.m.s. I — euj siii(j)SX 16X2/2?: N , both with and with- 
out TRS. In the opposite case of single-channel leads 
(N = 1), the distribution is highly non-Gaussian, see 
Fig. [| It has a logarithmic singularity (cusp) at zero 
current in the presence (absence) of time-reversal sym- 
metry. The asymptotics for large |/| are given by 



P(I) oc 



I/I- 9 / 4 TRS, 



no TRS. 



(12) 



The scattering matrix formula (g) allows us not only to 
find the statistical distribution of the pumped current I, 
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FIG. 3. (a) Distribution of current I — uij sin <j> SX1SX2 
for single-channel point contacts. The case of presence (ab- 
sence) of time-reversal symmetry is shown solid (dashed), (b) 
The same, but with capacitive interactions taken into account 
by a self-consistent Hartree approach, in the limit C <C pe 2 . 

but also the statistical correlation between / and the con- 
ductance G of a chaotic quantum dot. The correlation 
between / and G shows a remarkable dependence on the 
presence or absence of time- reversal symmetry (TRS): 
In Ref. 

it was shown that the statistical distribution 
of Rx x and Rx 2 is correlated with that of the conduc- 
tance in the presence of TRS only. Therefore, without 
TRS, / and G are statistically uncorrelated for a chaotic 
quantum dot, while they are correlated in the presence 
of TRS. In the latter case, we find that the width of the 
current distribution at a fixed conductance G is propor- 
tional to G 1 ' 2 for single-channel leads.tS 

Our main result, Eq. (||), is readily extended to in- 
clude the effect of a capacitive interaction in the email-, 
turn dot within a self-consistent Hartree treatment 0't£l 
Following Ref. |l^, the effect of the capacitive interaction 
is described by a self-consistent electric potential U. The 
potential U is related to the (kinetic) energy E and the 
Fermi energy Ep via Ep = E + U . Variation of X\ and 
X2 will cause a change of U and hence of E = Ep — U. 
(The Fermi energy Ep is kept constant.) Hence we have 
to deal with a simultaneous variation of E, X±, and X^. 
These simultaneous variations can still be described by 
the scattering matrix formula (§) provided we replace the 
response matrices Rxj {j = 1,2) by 



Rxi — > Rxi + Re 



dE 
dX~ 



Re 
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The derivative dE/dXj readal^ 
dE tvR Xl 



dXi 



irC/e 2 + trR E 



, i = 1,2 



(13) 



(14) 



where C is the geometrical capacitance of the dot. For 
many-channel contacts, inclusion of the interactions has 
no effect on the current distribution. For single-channel 
leads, we have computed the current distribution for the 
experimentally relevant case C <C e 2 p, where p is the 
(average) density of states in the dot, using the distri- 
bution of the matrices-ffxj and Re in the presence of 
capacitive interactions. cj The result, which is not much 



different from the non-interacting case, is shown in Fig. 
|b. 
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